Designing ferromagnetism in vanadium-oxide based superlattices 



Hung T. Dang-'^ and Andrew J. Millis^ 

'Department of Physics, Columbia University, 538 West 120th Street, New York, New York 10027, USA 

(Dated: June 4, 2013) 

Motivated by recent reports (Phys. Rev. B80, 241102) of room-temperature ferromagnetism in 
vanadium-oxide based superlattices, a single-site dynamical mean field study of the dependence of 
the paramagnetic-ferromagnetic phase boundary on superlattice geometry was performed. An exam- 
ination of variants of the experimentally determined crystal structure indicate that ferromagnetism 
is found only in a small and probably inaccessible region of the phase diagram. Design criteria for 
increasing the range over which ferromagnetism might exist are proposed. 

PACS numbers: 73.21.Cd,71.28.-f d,75.10.Lp 



I. INTRODUCTION 

"Materials by design" , the ability to design and create 
a material with specified correlated electron properties, is 
a long-standing goal of condensed matter physics. Super- 
lattices, in which one or more component is a transition 
metal oxide with a partially filled d-shell, are of great cur- 
rent interest in this regard because they offer the possi- 
bility of enhancing and controlling the correlated electron 
phenomena known^ to occur in bulk materials as well as 
the possibility of creating electronic phases not observed 
in bulk.'^ Following the pioneering work of Ohtomo and 
HwangjIS heterostructures and heterointerfaces of tran- 
sition metal oxides have been studied extensively. Ex- 
perimental findings include m etal -insulator transitions,* 
superconductivity,'^ magnetisnPEl and coexistence of fer- 
romagnetic and superconducting phases.'^ 
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FIG. 1: Ferromagnetic-paramagnetic phase diagram for 
La/SrVOs solid solution in plane of carrier concentration 
(changed by Sr concentration) and tilt angle in Pnma struc- 
ture but with all three Glazer's angles nearly equal. Dashed 
line indicates relation between carrier concentration and ro- 
tation amplitude in physically occurring bulk solid solution. 
From Ref . llOj 



In this paper we consider the possibility that ap- 
propriately designed superlattices might exhibit ferro- 
magnetism. Our work is partly motivated by a recent 
reporiP of room-temperature ferromagnetism in super- 
lattices composed of some number m of layers of LaVOa 
(LVO) separated by one layer of SrVOa (SVO), even 
though ferromagnetism is not found at any x in the bulk 
solid solution Lai_a;Sra;V03. Our study is based on a pre- 
vious analysiJ^ of the possibility of obtaining ferromag- 
netism in variants of the crystal structure of bulk solid 
solutions of the form Lai-ajSr^^VOa. A key result of the 
previous work was that ferromagnetism is favored by a 
combination of large octahedral rotations and large dop- 
ing away from the Mott insulating LaVOs composition. 
A schematic phase diagram is shown in Fig. [T] How- 
ever, as indicated by the dashed line in the figure, in the 
physical bulk solid solution, doping away from the Mott 
insulating concentration reduces the amplitude of the oc- 
tahedral rotations so that the physical materials remain 
far from the magnetic phase boundary. The motivating 
idea of this paper is that in the superlattice geometry, oc- 
tahedral rotation amplitude may be decoupled from car- 
rier concentration. The rotations can be controlled by 
choice of substrate while the carrier concentration can 
be controlled by choice of chemical composition and may 
vary from layer to layer of a superlattice. In effect, an 
appropriately designed superlattice could enable the ex- 
ploration of different paths in Fig. [ij 

In this study, we combine single-site dynamical mean 
field approximatior^ with realistic band structure cal- 
culations including the effects of the octahedral rotations 
to determine the ferromagnetic-paramagnetic phase dia- 
gram in supe rlattices with the crystal structures believed 
relevantfi^Ml ^^le experiments of Ref. El Unfortunately 
we find that the experimentally determined crystal struc- 
ture is in fact less favorable to ferromagnetism than the 
one found in the bulk solid solution, but we indicate 
structures that may be more favorable. 

The paper has following structure. The model and 
methods are described in Sec.|lll Sec. IIIII establishes the 
methods via a detailed analysis of the phase diagram of 
superlattices with no rotations or tilts. In Sec. |IV| we 
present the magnetic properties of superlattices with oc- 
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tahedral rotations similar to those observed experimen- 
tally. Section |V] is a summary and conclusion. 

II. MODEL AND METHODS 

A. Overview 

This paper builds on a previous study of the magnetic 
phase diagram of bulk vanadates. The new features rel- 
evant for the superlattices studied here are (i) the change 
in geometrical structure, including the differences from 
the bulk solid solution in the pattern of octahedral tilts 
and rotations and (ii) the variation of electronic density 
arising from superlattice structure. In the rest of this sec- 
tion we briefly summarize the basic theoretical method- 
ology (referring the reader to Ref . fTOl for details), define 
the crystal structures more precisely, explain the conse- 
quences for the electronic structure and explain how the 
variation of density appears in the formalism. 

B. Geometrical structure 

We study superlattices composed of layers of SrVOa 
(SVO) alternating with layers of LaVOa (LVO). If we ide- 
alize the structures as cubic perovskites, then the layers 
alternate along the [001] direction. In bulk, SVO crystal- 
lizes in the ideal cubic perovskite structure,'^^ while LVO 
crystallizes in a lower symmetry Pnma structure derived 
from the cubic perovskite via a four unit-cell pattern of 
octahedral tilts.^^ The crystal structure of bulk solid so- 
lutions Lai_a;Sra;V03 interpolates between that of the 
two end-members with the rotation amplitude decreas- 
ing as X increases. In the superlattice, the presence of a 
substrate and the breaking of translation symmetry can 
lead to different rotational distortions of the basic per- 
ovskite structure and also to a difference between lattice 
constants parallel and perpendicular to the growth direc- 
tion. 

Octahedral rotations in perovskites can be described 
using Glazer's notation. In the coordinate system de- 
fined by the three V-O bond directions of the original 
cubic perovskite, there are 3 tilt angles a, /3 and 7 with 
corresponding rotation axes [100], [010] and [001]. The 
tilt is in-phase if successive octahcdra rotate in the same 
direction, and anti-phase if they rotate in opposite di- 
rections. Rotational distortions of the cubic perovskite 
ABO3 structure may be denoted by a^Wc^ where fc 
can be — or denoting in-phase, anti-phase or no tilt- 
ing, respectively and a = a,b = f3,c = --y.LMMI Bulk LVO 
(Pnma) is of the type a^5+a^ with a — "/ — 8.7° and 
^ = 7.9° pEU 

For superlattices, substrate-induced strain may change 
the situation in a way w hich depends on the growth direc- 
tion. Experiment^^2llll confirm that the growth direction 
for the experimentally relevant superlattices is [001] (in 
the ideal cubic perovskite notation) and we focus on this 



case here. Recent experimental studies of superlatticeJl^l 
and of LaVOs thin films, which apparently have the same 
growth direc tion,!^ suggest that the rotations are of the 
type a-a+c~IiIliSI and indicate that the dominant rota- 
tion is around the axis defined by the growth direction: 
a ~ f3 3° and 7 w 11.5°. This distortion pattern is 
different from that occurring in bulk. To explore its ef- 
fects we set a = /3 = 3° and consider the consequences 
of varying 7. 

In bulk Lai_3;Srj:V03, while the 4-sublattice Pnma 
structure implies a difference in lattice constants, all V- 
O bond lengths are the same.^^ The difference in lat- 
tice constants arises from a difference in tilting pattern. 
Superlattices are typically grown on a substrate, and in 
epitaxial growth conditions the lattice constants perpen- 
dicular to the growth direction (which we denote here by 
a) are fixed by the substrate, while the lattice parameter 
along the growth direction (c) is free to relax. The result 
is a c/a ratio typically 7^ 1 contributed by both tilting 
and anisotropy in V-O bond lengths and possibly varying 
from layer to layer of the superlattice. F or th e experi- 
mentally studied superlattices, c/a ~ 1.02.'£'i2J The V-O 
bond lengths have not been determined but, as discussed 
in more detail in the Appendix, our studies indicate that 
all V-O bonds have essentially the same length. Further 
we show that a few percent differences have no signifi- 
cant effect on our study of ferromagnetism. In the rest 
of the paper we therefore ignore these distortions, setting 
all V-O bond lengths to be equal. 



C. Electronic structure 

We study superlattices designed to be similar to the 
system studied in Ref. |S1 In these superlattices, units of 
m layers of LaVOs are separated by one layer of SrVOa. 
To define the superlattice, we begin from LVO in the 
appropriate bulk structure, then break translation in- 
variance along the [001] (^-direction) by replacing every 
(m-f 1)**^ LaO plane with an SrO plane. Fig. [2^ shows 
such a superlattice with m = 3. 

We assume that the superlattice is grown epitaxially so 
that in-plane bond lengths and other aspects of the local 
structure including rotations are the same for all layers. 
We therefore take the electron transfer integrals which 
define the band structure to the be same for all layers. 
In this case the electronic structure of a superlattice is 
defined by adding the electrostatic potentials of the Sr 
and La ions to the basic translationally invariant hopping 
Hamiltonian describing the bulk materials. 

In our calculations we follow the common practice in 
studies of early transition metal oxides by assuming that 
the energy splitting between transition metal d-bands 
and oxygen p-bands is large enough to justify the use 
of a "frontier orbital" model focusing on the p-d anti- 
bonding bands which are mainly composed of vanadium 
t2(;-symmetry d-states. 
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FIG. 2; (Color online) (a) Schematic of superlattice lattice 
structure (LaV03)m(SrV03)i with m — 3. Vanadium sites 
indicated as circles with charge density indicated by shading: 
heavy shading (black online) indicating higher charge density 
and light shading (yellow online) indicating lower charge den- 
sity. LaO and SrO planes are shown as solid and dashed lines 
respectively. Nearest neighbor (t) and next-nearest neighbor 
{t') hoppings between vanadium sites indicated by arrows. 
The numbers on the right are VO2 layer indices, (b) Inset: 
pd-K hopping between t2g orbital and p-orbital. Main panel: 
two-dimensional nearest neighbor hopping t made of two pdn 
hoppings from xy orbital of one vanadium site to oxygen px 
or py orbital, then to xy orbital of another vanadium site. 



The Hamiltonian for the superlattice is thus 



H = Hkin + Ho 



H. 
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(1) 



where Hcouiomb describes the electron-ion interaction and 
electron-electron interaction between different sites and 
Honsite describes the d-d interactions, which we take 
to be on-site. Hum is a tight binding model, derived 
by using maximally-localized Wannier function (MLWF) 
technique^i^ to fit the t2g-derived antibonding bands. 
The detailed procedure is described in our previous 
work.liSl 

The kinetic Hamiltonian has the quadratic form 



k,Q,/3,CJ 
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where cj^^^ and ci^pa are electron creation and annihi- 
lation operators in reciprocal space with wavevector k. 
a and /3 are orbital and layer indices, and a is the spin 
index. 

We assume that the interaction takes the standard 
Slater-Kanamori forrrPSHU! which following Ref. [TU] we 
write as 



{U ~ 2 J) ^ riic,^nji3y 

+ iU-3J) J2 



(3) 



where the values of the on-site interaction U and the 
Hund's coupling J are U — 6eV ^ 22t and J = leV so 
that LVO is an insulator in bulk while SVO is a metal. 

In the approximation employed here, the superlattice 
is defined by the Coulomb interaction between the La/Sr 
ions and electrons. This, and the off-site part of the 
electron-electron interaction is containecPSl in 
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To construct Hei-iom we assume that the whole ion 
charge of SVO or LVO unit cell comes into the Sr or 
La site. Consider SrVOa, the valence of V is -1-4 (d^). 
If this one d-electron is removed, the SVO unit cell will 
have charge hence, in our model, Sr site has charge 
-1-1. Similarly, LaVOa has V+^ (d^), thus La site has 
charge +2. As a result, i?e/-ion has the form 



H, 



el — ion 
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where is electron-occupation operator at V-site z, e is 
the relative dielectric constant. The part Hf,i-ei is the 
inter-site Coulomb interaction of vanadium d-electrons 
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Hei-ei is treated in the Hartree approximation. Note 
that in Eq. (|6|, hi is the operator giving the total d- 
electron occupation of site i, while rij — (hj) is the ex- 
pectation value of d-electron occupancy at site j, which 
is determined self consistently. From H coulomb, the 
Coulomb potential Vi for site i is calculated using Ewald 
summation.^"* 

The dielectric constant e is an important parameter in 
Eqs. ([5j|6|. It accounts for screening on the scale of a lat- 
tice constant so bulk measurements are not directly rele- 
vant and an appropriate value has not been determined. 
Values ranging from 4 to 15 have been reported in the 
literature for similar systems .'^SISH Because the appropri- 
ate value of e has not been determined, we have studied 
several cases and present results mainly for e = 8, 15. 
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D. Methods 

We treat the on-site interaction terms using single-site 
dynamical mean field theory (DMFTjSl with the hy- 
bridization expansion continuous time quantum Monte 
Carlo (CTQMC) solveri^SThe superlattice effect is taken 
into account by the Coulomb potential V. We use the 
superlattice dynamic al me an field theory introduced by 
Potthoff and Noltin#S129| the form given in Ref. [H 
Here each V site i has a self energy (site local but de- 
pendent on site) determined from the solution of a quan- 
tum impurity model which has parameters fixed by the 
DMFT self-consistency equation linking the site local 
term of the lattice Green function {Ju to the quantum 
impurity model Green functior 



where 



. 47reeo i?i - Rj\ ^ 47reeo i?^ - Rsr\ 
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(8) 
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is a site dependent quantity, diagonal in spin and or- 
bital indices but linking different sites, derived from 
Eqs. ([5j[6|. The layers are coupled by a self-consistency 
condition which as discussed in Refs. I23I28I29 1 fixes 
both the hybridization function of the quantum impu- 
rity model and the layer-to-layer variation in the charge 
density. 

As described in Ref.fTUl it is advantageous to perform a 
site-local rotation to align the orbital basis to the local V- 
O bond directions of each octahedron before solving the 
impurity model. This reduces the sign problem in the 
CTQMC impurity solver and restores in-plane transla- 
tion invariance in the sense of making the self-consistency 
equations the same for all sites in a given plane. 

In a superlattice composed of N layers, it is in princi- 
ple necessary to solve N dynamical mean field problems, 
coupled by the self-consistency condition. However, we 
find (see section HI) that the susceptibility for a given 



layer of the superlattice may be determined from a bulk 
computation at the same local density and crystal struc- 
ture. Because the layer dependent density has no signif- 
icant dependence on the temperature or the many-body 
physics, it may be determined once from a band structure 
calculation and then bulk results with the appropriate 
density for a wide range of temperature may be used to 
infer the Curie temperature, substantially reducing the 
computational burden. 

The Curie temperature for ferromagnetism is deter- 
mined by extrapolating the inverse susceptibility x~^{J^^ 
to based on Curie- Weiss law ~ T — Tc. The test 
for the reliability of this method for has been done in 
Ref. Uni A similar approach can be found in literature.— 



III. RELATION BETWEEN SUPERLATTICE 
AND BULK SYSTEM CALCULATIONS 

In this section, we demonstrate that the magnetic 
phase diagrams of superlattice systems may be inferred, 
to reasonable accuracy, from the study of appropriately 
chosen bulk systems. This enables a considerable reduc- 
tion in the computation resources required. 
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FIG. 3: (Color online) Panel (a): Non-interacting density 
of states for bulk system at carrier density n = 1. Panels 
(b): Non-interacting density of states for different layers of 
(LVO)3(SVO)i superlattice for two different values of dielec- 
tric constant e = 5 (solid) and e = 15 (dashed) with hopping 
parameters t = 0.264eV and t' = 0.084eV. SrO plane is be- 
tween layers and 1 (the index is defined in Fig. [2|. The 
Fermi energy is at 0. 



We begin with a study of "untilted" or "cubic" su- 
perlattices: those in which all V-O-V bond angles are 
180°. We focus specifically on [001] superlattices in which 
the unit cell contains m layers LVO and one layer SVO, 
where to = 3, 4, 5. For orientation, we present the density 
of states (DOS) of the non-interacting system in Fig. [s] 
In obtaining these densities of states we used the simple 
tight binding parametrization. The DOS for the bulk sys- 
tem is shown in panel (a) . One sees the typical three-fold 
degenerate DOS for band, the Van Hove singularity 
is visible as a peak near the upper band edge. It is at 
high energy because the next-nearest neighbor hopping 
t' > 0. The remaining panels show the layer-resolved 
densities of states for the to = 3 superlattice. The upper 
two panels show layers sandwiched by La on both sides; 
the lower two panels show the layers adjacent to the SrO 
plane. The superlattice-induced changes in the density 
of states are seen to be relatively minor: the main effects 
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are a weak splitting of the van Hove peaks reflecting the 
breaking of translational invariance in the z-direction, 
and a relative shift in the positions of the van Hove peaks 
arising from band bending associated with the different 
charges of the Sr and La ions. 

Fig. [4] shows the layer-resolved charge density and in- 



verse susceptibility 



H 



plotted against temperature 



m{H) 

for three different superlattice structures corresponding 
to ni — 3,4,5. As expected from electrostatic consid- 
erations, the charge is lower for the VO2 planes nearer 
the SrO layer and the charge variation between layers is 
controlled by the dielectric constant. 

The magnetization m at the V sites on each layer was 
computed at field H = O.OleV / hb and the inverse sus- 
ceptibility was obtained as H/m. Linearity was veri- 
fied by repeating the computation using H = 0.02eV/^B 
(not shown). For the m = 3, e = 15 case (Fig. |4^), 
we extended the computation to the lower temperature 
T = 0.03eV; for the other two cases T = 0.06eV was 
the lowest temperature studied. The inverse susceptibil- 
ities are approximately linear in temperature at higher 
temperatures and in all cases, extrapolation to — 
reveals Tc < 0, implying absence of ferromagnetism. 

Especially for the layer nearest the SrO plane the 
curves exhibit weak upward curvature at the lowest tem- 
peratures studied. As shown in Ref. IIOI the curvature 
is a signature that the system is entering a Fermi-liquid 
coherence regime. The Fermi liquid coherence tempera- 
ture is highest for the layers nearest the SrO because the 
charge in these planes is farther from the n = 2 Mott 
insulating state. To verify this we followed Ref. [TO] and 
computed the Wilson ratio Rw for each layer of the su- 
perlattice for the case m = 3, e — 15, finding (not shown) 
that for each layer the Rw extrapolates to 2 at low tem- 
perature. The approach to the low temperature value 
is faster for layers with low density (near SrO planes) 
than for layers with high density (far from SrO planes). 
Rw = 2 is the value for a Kondo lattice, while ferromag- 
netism is characterized by an Ryi^ > 2.ESi We therefore 
believe that for "untilted" superlattices, the differences 
in x~^ among layers arise from differences in quasipar- 
ticle coherence scale, there is no evidence for ferromag- 
netism in this system, consistent with the solution of the 
corresponding bulk problem. 

To gain insight into the physics underlying the layer 
dependence of x~^ have computed x~^(T) for the 
cubic bulk system [Hcouiomb = 0, H}~in is constructed 
from the two-dimensional dispersion e(k) = —2t{cosk,j. + 
cos ky) — 4:t' cos kx cos ky) for carrier densities equal to 
those on the different VO2 layers. In Fig. |4j we present 
bulk calculations for n = 1.62 and n = 1.88 correspond- 
ing to the densities calculated for layer and 2 of the 
superlattice for all cases m = 3,4,5. For n — 1.88, bulk 
X~^ a,t T — 0.06, 0.10 and 0.14eV are very close to those 
of i = 2 layer of m = 3 superlattice, which has the 
same density. For m = 4, 5 superlattices, bulk n = 1.88, 
X~^{T) (not shown) almost coincides with those of i = 2 
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FIG. 4: (Color online) Temperature-dependent layer-resolved 
inverse magnetic susceptibilities for symmetry-inequivalent 
layers of untilted (LVO)m(SVO)i superlattice structures with 
different numbers of LVO layers m — 3,4 and 5. Layer is 
adjacent to SrO and layers 2 and 3 are between two LaO lay- 
ers. The relative dielectric constant is e = 15, magnetic field 
H = OMeV/tiB- The X~^iT) obtained from solution of bulk 
cubic systems with charge density set to the density on the 
given layer are also shown. "Bulk LO" ("BulkL2") denotes a 
calculation performed for a bulk system with density the same 
as for L = (L = 2) layer density. Inset: the electron layer 
density distribution corresponding to the susceptibility plot, 
a;-axis is the layer index, i/-axis is the layer density. On-site 
interactions U = 6eV, J = leV. 



layer. For bulk n — 1.62, the difference between bulk 
and superlattice L — layer is small. These calculations 
demonstrate a general rule: within the single-site DMFT 
approximation, the layer-resolved properties of a super- 
lattice correspond closely to those of the corresponding 
bulk system at a density equal to that of the superlattice. 

The superlattices of experimental relevance have crys- 
tal structures which are distortions of the "untilted" one, 
involving in particular a P2i/m structure characterized 
by a rotational distortion of the a~a~^c~ typ^l^^ involv- 
ing a large rotation about an axis approximately paral- 
lel to the growth direction and much smaller rotations 
about the two perpendicular axes. Fig. [5] compares the 
non-interacting DOS of bulk and (LV0)3(SV0)i super- 
lattice systems (both with the same P2i/ni structure) 
calculated using DFT and a MLWF parametrization of 
the frontier bands. The DOS of bulk system is shifted so 
that it has the same carrier density as layers of the su- 
perlattices near SrO plane. For three different structures 
(untilted structure and P2i/m structure with 7 = 11.5° 
and 16°), the basic features of the partial DOS are simi- 
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FIG. 5: (Color online) Comparison between bulk LVO partial 
DOS (positive curves) and (LaV03)3(SrV03)i superlattice 
layer DOS of layers near SrO (negative curves) derived from 
band structure calculations (DFT-fMLWF). Both systems 
have the same lattice structure for each case: untilted struc- 
ture for the top panel and P2i jm structure (Glazer's notation 
a~a"'"c~) with q = /3 = a = 3° and 7 = c = 11.5° and 16° 
for other panels. The DOS of bulk system is shifted towards 
higher energy so that bulk carrier density is the same as layer 
density of superlattices for the layers near SrO (n ~ 1.55). 
The vertical dashed line marks the Fermi level. 
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FIG. 6: (Color online) Comparison in temperature depen- 
dent inverse susceptibility between bulk LVO (solid lines) and 
(LaV03)3(SrV03)i superlattice (dashed lines). Both have the 
same lattice structure P2\/m with tilt angle 7 = 11.5, 16 and 
18°. Bulk system has the same densities as those of layers 
of superlattice near and far from SrO planes (n = 1.55, 1.95). 
Left column: the plots in wide temperature range. Right col- 
umn: the expanded views near zero temperature. 



lar between bulk and superlattice. The translation sym- 
metry breaking in z-direction leads to small extra peaks 
in the superlattice DOS. These differences are smoothed 
out by the large imaginary part of the DMFT self en- 
ergy. Because the DMFT equations depend only on the 
density of states it is reasonable to expect that, as in the 
untilted case, they will therefore give the same results in 
the superlattice as in the bulk material with correspond- 
ing density of states. 

To verify that this is the case we have also compared 
bulk and superlattice susceptibilities for tilted structures. 
The four VOg octahedra in a unit cell are related by ro- 
tation, so an appropriate choice of local basis means that 
only one calculation needs to be carried out for a given 
layer. Fig. [6] compares the inverse susceptibilities for an 
TO = 3 superlattice to calculations performed on a bulk 
system with the same P2i/to structure. In these calcu- 
lations, we choose 7 = 11.5°, 16°, 18° and dielectric con- 
stant e = 8. We see that in this case, as in the "untilted" 
case, the superlattice inverse susceptibilities x~^{T) are 
almost the same as those for bulk system calculated at 
the same density, with differences only resolvable in the 
expanded view for the largest tilt angles. 



IV. SUPERLATTICES WITH GdFe03-TYPE 
ROTATION 

In this section we present and explain our results for 
the magnetic phase diagram of (LVO)„i(SVO)i super- 
lattices with the P2i/to structure (Glazer's notati on 
a~a+c~) reported for the experimental systems .^i^Ell In 
these structures in-plane rotation along the growth di- 
rection z = [001] is large 7 = 0= 11.5° (presumably 
because of the strain imposed by the substrate), while 
the out-of-plane rotation is small [a — (3 — a — per- 
haps because the system is free to relax along the growth 
direction. We concentrate on the effect of the large rota- 
tion by fixing the in-plane angles to 3° while varying the 
out-of-plane angles over a wide range from 10° — )■ 18°. 

Based on the results of Section |III| we generate a phase 
diagram for the superlattice from calculations for a bulk 
system which is a P2i/ to distortion of the ideal cubic per- 
ovskite structure of chemical composition LaVOa. The 
bulk system results are presented as a phase diagram in 
the plane of carrier concentration and 7-rotation. Spe- 
cific layers of the superlattice will correspond to particu- 
lar points on the phase diagram, with the layer dependent 
density fixed by number of LVO layers to and the dielec- 
tric constant e and the rotation fixed by the substrate 
lattice parameter. 



7 




-0.5 0.0 0.5 
energy (eV) 



1.5 



FIG. 7: (Color online) Partial DOS derived from 
DFT+MLWF for "bulk" P2i/m structure (Glazer's notation 
a~a^c~) with a = P = a = 3° and 7 = c changing from 
10° to 16°. Only fog bands are plotted because eg bands are 
negligible in this range of energy. 



We use DFT+MLWF methods to obtain the frontier 
orbital band structure for the f2g-derived antibonding 
bands. Fig. [7] presents representative results for the or- 
bitally resolved local density of states. In this figure the 
orbitals are defined with respect to the local basis defined 
by the 3 V-0 bonds of a given VOg octahedron. We de- 
fine z = [001] as the axis (approximately parallel to the 
growth direction) about which the large rotation occurs. 
Fig. [7] shows that yz and zx orbitals are almost degen- 
erate, while xy orbital is strikingly different. The DOS 
of xy orbital maintains the shape of a two-dimensional 
energy dispersion with a van Hove peak well above the 
chemical potential, similar to the bulk cubic structure 
(see e.g. Fig. [3^). There are noticeable differences only 
at very high rotation angles. On the other hand, yz and 
zx orbitals are spread out with two small peaks, because 
hoppings along x or y directions (more distorted) are 
different from those along z-direction (less distorted). 
When the distortion gets larger, the d-bandwidth be- 
comes smaller, the xy peak gets larger and slightly closer 
to the Fermi level, and yz and zx peaks near the Fermi 
level also develop. 

Based on Hband{^) generated by DFT+MLWF, we 
carry out DMFT calculations for in-plane rotation an- 
gle 7 to get curves whose extrapolations define the 
Curie temperatures Tc- Fig.[8]shows how Tc evolves when 
the rotation angle 7 increases from 10 to 18°. In this fig- 
ure, we consider two different carrier densities n = 1.55 
and 1.95, corresponding to the band structure prediction 
for the layer densities of layers near and far from SrO 
planes in the superlattice. Tc for n = 1.95 is a slow func- 
tion of rotation and is always negative for the range of 
7 under consideration, while Tc for n = 1.55 increases 



faster, so that the system becomes ferromagnetic when 
7 is between 14 and 16°. Ferromagnetism is therefore 
expected only in superlattices with very large rotations, 
and then only in the layers with large hole doping (i.e. 
the layers closest to the SrO planes). 
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FIG. 8: (Color online) Inverse susceptibility x"'^ vs. tem- 
perature T for bulk P2^/m structure of LaVOs at densities 
n — 1.55 (black circle solid lines) and n = 1.95 (red diamond 
dashed lines) for rotation angle 7 increasing from 10 — >■ 18°. 
On-site interaction U = 6eV and J = leV. Left column: the 
circles and diamonds are data points, the solid and dashed 
lines are fitted from these data points. Right column: ex- 
panded view at small region. The vertical dashed line 
marks zero temperature. 

From a range of calculations such as those shown in 
Fig. [8] we have constructed the superlattice magnetic 
phase diagram shown in Fig. |9] Similar to Ref. [lOl there 
are uncertainties in our extrapolation for Curie tempera- 
ture, we consider 0.004eV as the error bar for positions on 
the phase diagram. Thus, Tc < 0.004eV is considered as 
Tc = within the error bar. We see that ferromagnetism 
is favored only for very large rotations, much larger than 
the 11° determined experimentally, and only for carrier 
concentrations far removed from n — 2. We may locate 
the experimentally studied superlattices on this phase 
diagram. For an m = 3 superlattice, band structure cal- 
culations indicate layer densities 1.55 for layers near SrO 
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plane and 1.95 for the other layers. The experimentally 
determined rotation angle is ~ 11.5°. These two points 
are indicated by squares in Fig. [9] 
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FIG. 9: (Color online) The magnetic phase diagram with x- 
axis carrier density n and j/-axis tilt and rotation angle along 
z — [001] direction 7 for bulk system LVO with the same 
type of distortion as for (LVO)m(SVO)i superlattices {P2\/m 
structure), in-plane tilt angles a,/3 « 3°. On-site interactions 
U — 6eV, J = leV. The white regime indicates absence of 
ferromagnetism (Tc < 0.004ey), the colored regime indicates 
ferromagnetism with Tc indicated by the color bar. Also in- 
dicated are results for bulk Lai_2,Sri,V03 in the Pnma struc- 
ture, from Ref. llOl Note that in the calculations for the Pnma 
structure all three tilt angles are almost the same. 

It is interesting to compare our results to those previ- 
ously obtainecf^for the bulk solid solution Lai-j^Sr^jVOs 
{Pnma structure). The dashed line in Fig. [O] shows the 
theoretically estimated phase diagram for the bulk solid 
solution. We see that the bulk structure is more favorable 
for ferromagnetism than the superlattice structure. An 
important difference between the Pnma structure and 
the P2i/m of the superlattice is that in the former case 
all three tilt angles are of comparable magnitude whereas 
in the P2i/m structure only one rotation is large. We be- 
lieve that this difference is responsible for the difference 
in phase boundary. 



V. CONCLUSIONS 

In this paper, we have studied the possibility of ferro- 
magnetism in superlattice structures of vanadium oxides 
derived from LaVOa and SrVOs. Our investigation was 
based on the idea that ferromagnetism depends on an in- 
terplay between carrier density and octahedral rotation, 
and while these are coupled in bulk (see the solid solution 
curve in Fig. |9]) they may be decoupled in the superlat- 
tice. In particular, the charge density varies across the 
superlattice, being lowest near the SrO planes, while the 



rotation angle is controlled by the substrate. Thus in 
an appropriately designed superlattice at least some por- 
tions of the system might be moved closer to (or perhaps 
into) the ferromagnetic region. In several important as- 
pects this idea is consistent with calculations. We find 
that the local carrier density determines the local mag- 
netic susceptibility (see section III I and the density/tilt 



angle relationship may be significantly altered (see solid 
line and square points in Fig. [9]). 

However, we find that the P2i/m octahedral rota- 
tion pattern characteristic of experimentally discovered 
superlattices is in fact less favorable to ferromagnetism 
than the Pnma pattern characteristic of bulk materials 
(compare the phase boundaries in Fig. [9]). Thus while 
the general idea that an appropriately designed super- 
lattice might provide conditions favorable for ferromag- 
netism thereby providing a potential explanation for the 
remarkable experimental report of room-temperature fer- 
romagnetism in (LaV03)m(SrV03)i superlattices with 
TO = 3, 4, 5, 6 by Liiders et. al.,'^ (even though there is no 
ferromagnetism in the bulk solid solution), our detailed 
findings arc not consistent with the experimental result. 

Our results indicate that designing ferromagnetism 
into a vanadate superlattice will require both large ampli- 
tude rotations about the growth axis and also substantial 
rotations about the other two axes. Rotations about the 
growth axis arise from substrate-induced strain, so choos- 
ing substrates with smaller lattice parameter would be 
desirable. Introduction of rotations about the orthogo- 
nal axes may be done by replacing the La with a smaller 
counterion such as Y. 

Our study has certain limitations. The calculations 
employ a frontier orbital model which includes only the 
i2g-derived antibonding bands. DFT-I-DMFT calcula- 
tions based on correlated atomic-like d-states embedded 
in the manifold of non-correlated oxygen states provide 
a more fundamental description. Our previous work?^ 
indicates that the two models give very similar results 
if both calculations are tuned so that bulk LaVOs is 
a Mott insulator, but the implications of the full (but 
computationally very heavy) DFT-I-DMFT procedure for 
the superlattice problem remain an open problem for fu- 
ture research. Further, our calculations are based on the 
single-site DMFT approximation, which includes all local 
effects but misses inter-site correlations. While it is gen- 
erally accepted that these calculations give the correct 
trends and qualitative behavior, the quantitative accu- 
racy of the methods is not known. Unfortunately, as yet 
cluster extensions of DMFT are prohibitively expensive 
for the multiband models considered here. 

The experimental results of Liiders et. alW there- 
fore provide an interesting challenge to materials the- 
ory. They indicate that superlattices display ferromag- 
netism when the corresponding bulk solid solutions do 
not, whereas the present state of the art of real materials 
dynamical mean field calculations suggests that super- 
lattices should be less likely to display magnetism than 
the corresponding bulk solid solutions. This discrepancy 
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requires further investigation. 



bond angles are compatible with the experiment. 
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Appendix: Lattice constant and V-O bond length 
ratio 

In this appendix, we present a more complete discus- 
sion of the strain-induced lattice distortions. The in- 
plane lattice constant of a superlattice epitaxially grown 
on a substrate matches that of the substrate and may 
therefore be different from the lattice constant preferred 
in a free-standing film or bulk material. The out-of-plane 
lattice constant is typically free to relax, and in the pres- 
ence of an in-plane strain may also be different from that 
found in bulk materials. 

A difference in V-V distance may arise from a change 
in V-O bond length or from a difference in buckling of 
V-O bonds. We consider both possibilities here, but first 
remark that the main differences in structure between 
bulk and experimentally studied superlattices arise from 
differences in octahedral rotation. In the experimentally- 
studied superlattices, the in-plane V-V distance is in fact 
slightly less than the V-V distance in LVO. The V-O bond 
lengths have not been measured for the superlattice, but 
to a high degree of accuracy we are able to reconstruct 
the measured superlattice using the measured tilt an- 
gles given from experimentJ^^ structure, assuming that 
all V-O bond lengths are equal. Assuming the P2i/m 
structure, we varied the in-plane and out-of-plane V-O 
bond lengths to fit the experimental data and found that 
c/a 1.02 only when the mean bond length d is found 
in the range from 1.983 to 2A depending on which ex- 
perimental result is fit but in all cases the V-O bond 
lengths are found to be equal to within an accuracy 
of 0.3%. Therefore, we believe that all the V-O bond 
lengths should, to a good approximation, be the same. 
The structure used in our calculations is presented in Ta- 
ble, [l] Although there are slight mismatches in in-plane 
angle and lattice constants, the c/a sa 1.02 ratio and 



1.5 
1.0 
0.5 
0.0 



-1 ' ' 1 

Ay \y£.) 


'l 
'l 
'l 

lit 

\ 1 
1 1 

I \ 1 


(a) 






(b) 


,v 

IS 

1 1 

' \ A 



-1.0 



-0.5 0.0 0.5 1.0 
energy (eV) 



1.5 



FIG. 10: (Color online) Partial DOS for bulk LVO with c/a = 
1.02 with the c/a ratio due to a change in V-O bonds (panel 
(a)) and to P2\/m lattice structure with a = /3 = 0, 7 = 11.5° 
(similar to superlattice structure) (panel (b)). The dashed 
blue curve is the xy orbital, the solid red curve is the de- 
generate yz (or zx) orbital. The dashed vertical line marks 
the Fermi level. 

Changing the amount of rotation has a different effect 
on the electronic structure than does changing the ra- 
tio of V-O bond lengths. Fig. 10 compares the partial 



DOS for the two cases, using as example a hypotheti- 
cal LaVOa crystal with c/a — 1.02. The upper panel 
presents the DOS for the untilted structure with straight 
V-O-V bonds and the c/a ratio induced by a diflFerence in 
in-plane and out-of-plane V-O bond lengths. The lower 
panel presents the case of all equal V-O bonds, with the 
c/a ratio produced by octahedral rotations about the z 
axis. The densities of states are quite different, but can 
be understood from the simple energy dispersion 



rj/(k) 



2t\\ (2 — cos kx — cos ky)- 

+ 4i'| (1 — cos fcj; cos fcy), 

(k) = 2i||(l - cosfc^)+ 

+ 2i_L(l - cos k^) +4t^(l 



cos kx cos kz), 



(A.l) 

where ty and t_L are the in-plane and out-of-plane near- 
est neighbor hopping integrals and tj| ^ are the second 
neighbor hoppings. The lower band edge is assumed to 
be the same for all orbitals but we assume that the lat- 
tice distortions lead to different values for the in-plane 
and out of plane hoppings. 

The lower band edge is defined to be zero and is inde- 
pendent of the distortion. The energy of the upper edge 
of the xy band is e^y (tt, tt) = Siy and of the xz/yz bands 
is ea;z(7r,7r) = 4t|| -I- At^. The positions of the van Hove 
singularities are at k = (0, tt) or (tt, 0). For xy band there 
is only one van Hove peak, at = 4<|| -I- %t'\ while for 
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TABLE I: WyckofF positions and lattice constants for LaVOs with P2\/m structure from our calculation based on tilt angles 
taken from the experiment (Ref. lT3|l . La positions are fixed manually but do not affect lattice constants or c/a ratio. In 
our notation, d is the V-O bond length, ao,&o and co are lattice constants, /3o is the angle between ao and co, a and c are 
pseudocubic lattice constants (growth direction is along c direction, a is perpendicular to c). 
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FIG. 11: (Color online) Inverse susceptibility vs. tempera- 
ture for cubic structure of bulk hole-doped LVO. The in-plane 
and out-of-plane bondlengths are changed so that the octahe- 
dral volume is unchanged: tensile strain (c/a = 0.98 - black 
lines), no strain (c/a = 1.00 - red lines) and compressive strain 
(c/a = 1.02 - blue lines). Two levels of hole doping are con- 
sidered: n = 1.55 (solid lines) and n = 1.95 (dashed lines). 
These lines are linear fits for the data points. 



xz band, there are two van Hove peaks at — 



4t_L -|- 8^'. When <_l is different from t\\ 



and 
the differ- 



peaks of zx band e]^ — . 

With these definitions, we are in a position to under- 
stand the changes in the band structure. When the V- 
O bond lengths change (Fig. 10 1) so that the z-bond is 
longer and the in-plane bond is shorter but the octahe- 
dral volume is unchanged, the band structure calculation 
indicates that t_L decreases but increases slightly. The 
difference between the bandwidth of the xy and xz band- 
widths is 4i|| — 4tj^ which is the same as the splitting 
between the van Hove peaks in the xz/yz bands. On 
the other hand, if the c/a ratio is produced by rotation, 
(Fig. 10 d), the change is opposite. The in-plane hop- 



ping tj| decreases because of the buckled in-plane V-O-V 
bonds, while the out-of-plane hopping t±^ is unchanged. 
The xy band therefore narrows substantially relative to 
the xz/yz bands. In addition the splitting of the van 
Hove peaks is greater. From the bandwidth of xy and 
zx bands (Fig.fTolD), iy « 0.225eV, i_L ~ 0.35eV, the van 
Hove peak distance is sa 0.5eV, which is compatible with 
the peak positions shown in Fig. |10[ 3. 

We tested with DMFT calculations for the Curie tem- 
peratures with the V-O bondlength changed. Fig. 11 is 



ence in bandwidth of xy and zx orbitals is 4t|| — 4ij^, 
which is also the distance between the two van Hove 



the temperature-dependent inverse susceptibility derived 
from DMFT for the bulk cubic structure with the c/a 
ratio changing from 0.98 (tensile strain) to 1.02 (com- 
pressive strain). For all the levels of hole doping under 
consideration, the results are nearly the same for every 
case of c/a ratio. We conclude that even when the V-O 
bondlength changes within the physical range, the ferro- 
magnetism is not affected. However, we also found that 
when the V-O bondlength is such that c/a > 1.06 or 
< 0.90, there is large orbital polarization and the fer- 
romagnetism can be largely affected. But that range is 
unphysical and can be neglected in the context of this 
work. 
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